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We study the post-Newtonian expansion of a class of Lorentz-violating gravity theories that reduce to khrono- 
metric theory (i.e. the infrared limit of Horava gravity) in high-acceleration regimes, and reproduce the phe¬ 
nomenology of the modified Newtonian dynamics (MOND) in the low-acceleration, non-relativistic limit. Like 
in khronometric theory, Lorentz symmetry is violated in these theories by introducing a dynamical scalar field 
(the “khronon”) whose gradient is enforced to be timelike. As a result, hypersurfaces of constant khronon de¬ 
fine a preferred foliation of the spacetime, and the khronon can be thought of as a physical absolute time. The 
MOND phenomenology arises as a result of the presence, in the action, of terms depending on the acceleration 
of the congruence orthogonal to the preferred foliation. We find that if the theory is forced to reduce exactly 
to General Relativity (rather than to khronometric theory) in the high-acceleration regime, the post-Newtonian 
expansion breaks down at low accelerations, and the theory becomes strongly coupled. Nevertheless, we iden¬ 
tify a sizeable region of the parameter space where the post-Newtonian expansion remains perturbative for all 
accelerations, and the theory passes both solar-system and pulsar gravity tests, besides producing a MOND 
phenomenology for the rotation curves of galaxies. We illustrate this explicitly with a toy model of a system 
containing only baryonic matter but no Dark Matter. 

PACS numbers: 95.35.+d,04.50.Kd 


I. INTRODUCTION 

The year 2015 marks the hundredth anniversary of General 
Relativity (GR). This elegant theory has been greatly success¬ 
ful at interpreting and predicting gravitational phenomena on 
a huge range of length-scales, velocities, gravitational-field 
strengths and space-time curvatures. On small length-scales, 
submillimeter experiments verified the validity of Newto¬ 
nian gravity, to which GR reduces in the quasi-static weak- 
field regime characterizing these experiments, down to micro¬ 
meter scales HI El- Newtonian gravity has been historically 
tested in the solar system, but in the course of the twentieth 
century technological progress made it possible to test also the 
first post-Newtonian (1PN) corrections to Newtonian dynam¬ 
ics 011. i.e. the GR corrections of fractional order 0(v/c) 2 , 
with v being the system’s characteristic velocity. Indeed, 
these 1PN solar-system tests date back to the first triumph of 
GR, i.e. Einstein’s prediction of the correct perihelion shift 
for Mercury, and later came to include also light-deflection 
measurements, time-delay and gyroscopical-precession ex¬ 
periments, as well as exquisite verifications of one of GR’s 
building blocks, i.e. the equivalence principle. However, be¬ 
cause velocities in the solar system are v < 10 4 c, and the 
gravitational fields are weak (i.e. 4>n/c 2 < 10~ 6 , (f>N be¬ 
ing the Newtonian potential), tests of the GR dynamics be¬ 
yond this weak-field, mildly relativistic regime are impossible 
there. 

A glimpse at the workings of gravitation in a different 
regime is offered by binary pulsars, i.e. systems comprising 
of a pulsar (which allows accurate tracking of the orbital pe¬ 
riod), and another compact star (typically a neutron star or 


a white dwarf). These systems, the first of which was dis¬ 
covered in 1974 0. have velocities that are not much larger 
than in the solar-system (v < 10~ 3 c), but present large grav¬ 
itational fields/curvatures inside the compact stars. In this 
mildly relativistic but strong-field regime, GR predicts that 
gravitational waves (GWs) should be copiously emitted, thus 
carrying enough energy and angular momentum away from 
the binary to produce an observable backreaction on its or¬ 
bital evolution. Indeed, as the binary shrinks as a result of 
GW emission, its period should decrease. This effect has in¬ 
deed been observed in binary pulsars, and the period’s rate of 
change matches perfectly the GR prediction, thus providing 
indirect evidence of the existence of GWs EE). 

Finally, “advanced” ground-based GW interferometers, 
such as Advanced LIGO, Advanced Virgo and KAGRA, will 
come online in the next few years, and are expected to detect 
GWs directly before the end of this decade. Because the main 
GW sources for these detectors are expected to be binaries 
of neutron stars and/or black holes at small separations (and 
thus with relative velocities v ~ c), these interferometers will 
provide the first test of GR in the currently unexplored highly 
relativistic and strong field regime (see e.g. Refs. EUl for 
two recent reviews). 

Despite GR’s past triumphs and the busy experimental ac¬ 
tivity to test it even further with GWs, signs that something 
could be wrong with our understanding of gravity might al¬ 
ready be hidden in plain sight in cosmological data. In the last 
two decades, observations of the Cosmic Microwave Back¬ 
ground (CMB), type-la supernovae and the large-scale struc¬ 
ture of the universe pointed to the existence of a Dark Matter 
component and a cosmological constant (or a dynamical Dark 
Energy component); see e.g. Ref. ifTOl for a review. While this 
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“concordance” ACDM model is in agreement with essentially 
all observations so far, it is theoretically unappealing because 
“naturalness” arguments can explain neither the small value 
of the cosmological constant compared to the Planck scale, 
nor why it has only recently started to drive the expansion 
of the universe IfTTirnn . In the light of the ACDM model’s 
“unattractiveness”, it makes sense at least to ask the question 
of whether the existence of a Dark Sector may simply be an 
artifact of our use of GR to explain cosmological observa¬ 
tions. Because these observations are well within the weak- 
field, mildly relativistic regime tested in the solar system, the 
answer to this question would seem to be negative. This rea¬ 
soning, however, neglects some important considerations. 

First, the Newtonian and PN dynamics that are verified in 
the solar system and in binary pulsars are expansions around 
the Minkowski geometry. This is not suitable for describ¬ 
ing cosmological scales, which are rather described by the 
Robertson-Walker geometry (and by perturbative expansions 
around it). While in GR perturbative expansions around the 
two space-times behave in similar ways, the same is not guar¬ 
anteed to happen in more general gravity theories. For in¬ 
stance, certain gravity theories may have a screening mech¬ 
anism built in, which triggers modifications away from the 
GR behavior only under certain conditions mm, e.g. on 
large cosmological scales. It is remarkable that hints in favor 
of such a screening mechanism might be hidden in already 
available cosmological data. Indeed, observations of veloci¬ 
ties on galactic and galaxy-cluster scales seem to point at the 
existence of a universal acceleration scale a o = 1.2 x 10 -10 
m/s 1 2 ~ c/7 ( ) (where H 0 is the present Hubble rate). 

The appearance of such a universal scale is not an obvious 
feature of the ACDM model, which in order to interpret these 
data has to be supplemented with hypotheses about the bary- 
onic physics and its feedback on the growth of structures (see 
e.g. Refs. MM for recent reviews about galaxy formation 
in the ACDM model). Even worse, these additional assump¬ 
tions need to be finely tuned to correctly reproduce the data, 
at least in specific cases HMD- The appearance of a uni¬ 
versal scale linked to the Hubble rate fits instead in the logic 
presented above, in which deviations from the GR behavior 
appear when one moves away from perturbative expansions 
over Minkowski space toward expansions over a Robertson- 
Walker space-timePl Alternatively, one can devise gravity the¬ 
ories that include an acceleration-based screening mechanism, 
whereby GR is recovered in high-acceleration regimes (i.e. 
in the solar system and binary pulsars) and modified in low- 
acceleration ones, where the ACDM postulates the existence 
of Dark Matter (and Dark Energy). Indeed, the appearance 
of the universal acceleration a o in observations of galaxies 
and galaxy clusters may be a guiding principle in construct¬ 
ing a theory of gravity alternative to Newtonian theory/GR, 
in the same way in which Kepler’s laws were instrumental 
in overcoming the Aristotelian/Ptolemaic mechanics. These 
acceleration-based attempts, which are known under the name 


1 The Robertson-Walker geometry globally reduces to the Minkowski one 

when the Hubble expansion rate is zero at all times. 


of “Modified Newtonian Dynamics” (MOND) lf22W24ll . are 
not yet completely successful, because to explain observa¬ 
tions of galaxy clusters they still need some residual “dark 
missing baryons”, with mass roughly twice that of observed 
baryons lfl9l and possibly in the form of molecular hydro¬ 
gen l25l . [Note that this is not in contrast with the estimate 
of the baryon density coming from Big-Bang Nucleosynthe¬ 
sis (BBN), since about 30% of the baryons produced during 
BBN are still undetected, and only 4% are observed in clus¬ 
ters 1261 .1 Nevertheless, the appearance of a universal scale 
in the data is a genuine empirical feature, the explanation of 
which is still poorly understood. 

Another independent motivation for considering possible 
modifications of GR comes from its intrinsic incompatibil¬ 
ity with quantum field theory, i.e. the long-known fact that 
GR, when quantized, is not power-counting renormalizable 
in the ultraviolet (UV) regime, where it should be replaced 
by a (yet unknown) quantum theory of gravity. In addition, 
GR generically predicts the existence of curvature singulari¬ 
ties in time evolutions starting from regular initial data. Even 
though these singularities are conjectured to be always en¬ 
closed by black-hole horizons and thus inaccessible to out¬ 
side observers |[27l[28l . their existence is a disturbing feature 
that one expects should be solved by a full quantum theory of 
gravity. 

A candidate quantum-gravity theory that addresses these 
two problems is given by Horava gravity 12911301 . This theory 
breaks boost-symmetry (and thus Lorentz invariance) in the 
gravitational sector by adding to the action terms that are of 
fourth and sixth order in the spatial derivatives of the met¬ 
ric. In simpler scalar toy models, these terms are enough 
to achieve UV power counting renormalizability EH ED, 
and the hope is that the same will happen for spin-2 gravi¬ 
tons. Also, the presence of the higher-order terms in the 
spatial derivatives is expected to smooth the curvature sin¬ 
gularities typically forming in GR evolutions ll32l . On as- 
trophysical scales, Horava gravity is practically indistinguish¬ 
able from its low-energy limit, sometimes called “khronomet- 
ric theory” 1 301 l33l . This theory has been extensively stud¬ 
ied, thanks also to the fact that it is closely related P4ll35l to 
another previously introduced and actively scrutinized fam¬ 
ily of phenomenological boost-violating gravity theories, i.e. 
Einstein-TEther theories mm. Remarkably, khronomet- 
ric theory (and thus Horava gravity) has been shown to pass 
all experimental tests, i.e. submillimeter tests ||32l , absence 
of gravitational Cerenkov radiation l38l . solar-system ex¬ 
periments E3E3ED, binary- and isolated-pulsar observa¬ 
tions ED S3, and existence of regular black holes forming 
from gravitational collapse (so as to agree with astrophysical 
observations of black-hole candidates) 143)4491 , in regions of 
parameter space where khronometric theory is stable at both 
the classical and quantum levels. 

An attempt at modifying khronometric theory and Horava 
gravity to account for the presence of a universal accel¬ 
eration scale in galaxy and galaxy-cluster data was done 
in Ref. [50), which introduced a theory that reduces to a 
(very special) khronometric/Horava-gravity theory in high- 
acceleration regimes, and which produces a MOND behav- 
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ior in the low-acceleration, non-relativistic/weak-field regime 
relevant for galaxies and clusters. This theory clearly shares 
both the flaws and the blessings of MOND that we mentioned 
above, namely it accounts for the appearance of a universal 
acceleration without finely tuned baryonic physics/feedback, 
but may still need some form of Dark Matter in the center of 
galaxy clusters]^] Also, the theory of Ref. |50l is related to 
some of the older theories proposed to obtain a MOND-like 
phenomenology in the non-relativistic limit - namely tensor- 
vector-scalar gravity (TeVeS) ED and generalized Einstein- 
/Ether theories (5211531; c.f. also Ref. 01 for an extensive 
review of the theories giving a MOND phenomenology but 
is better motivated theoretically, because it reduces to a viable 
quantum gravity model such as Hofava gravity at high accel¬ 
erations. 

In this paper, we will work out the 1PN expansion of the 
theory of Ref. Il50l . in both the high- and low-acceleration 
regimes. We will show that if one imposes that the the¬ 
ory reduces to GR in the high-acceleration regime, a strong¬ 
coupling problem arises in the low-acceleration regime when 
1PN terms are considered in the dynamics, and this would ruin 
the agreement with the observed rotation curves of galaxies. 
Indeed, we will show that while these observations are repro¬ 
duced in the Newtonian limit, the 1PN dynamics is strongly 
coupled, as a result of which the 1PN terms become domi¬ 
nant over the Newtonian ones in regimes accessible by galaxy 
rotation curves. However, we will then show that a simple 
slight generalization of the theory of Ref. ESI allows us to 
avoid this strong-coupling problem, i.e. one can obtain a fully 
viable theory by relaxing the assumption that the dynamics 
should reduce exactly to GR in the high-acceleration limit. 
We will therefore end up with a theory that (i) presents a well- 
behaved (i.e. perturbative) PN expansion at all accelerations; 
(ii) passes submillimeter, pulsar and solar-system tests; (Hi) 
reduces to a general khronometric theory (and thus to Hofava 
gravity) at high accelerations; (iv) gives a MOND-like phe¬ 
nomenology at the low accelerations characterizing galaxies 
and clusters. 

This paper is organized as follows. In section [II] we intro¬ 
duce the theories under investigation. The dynamics of these 
theories in the high-acceleration regime, as well as the exper¬ 
imental/theoretical constraints on it, are discussed in section 


III The low-acceleration regime, and in particular the 1PN 


dynamics, is discussed in section [TV] both in the general case 
and for the special case of a galaxy accreting gas. We show 
that the low-acceleration 1PN dynamics is strongly coupled 
in a certain region of parameter space, and that this may jeop¬ 
ardize the agreement of the theory with data on the scales of 
galaxies. In section [V] we identify this region, and show that 
the theories that we consider remain viable in large portions 


2 The necessary amount of Dark Matter is smaller than in the ACDM model. 
Indeed, as mentioned above, it might be sufficient to identify this Dark Mat¬ 
ter with some of the “missing dark baryons” that are predicted by BBN, but 
which are not observed in the local universe in the form of visible matter. It 
has been proposed that these dark baryons may be in the form of molecular 
hydrogen 1251 . 


of the parameter space. A final discussion is then presented in 
section m 

We will also use a metric signature (—++), and we will 
denote space-time indices by Greek letters and spatial ones 
by Latin letters. Spatial vectors are also denoted by an over¬ 
arrow. We will set c = 1 throughout this paper, except when 
dealing with PN expansions in sections IV B||IV C and in the 
Appendix, where we reinstate the factors 1/c as PN book¬ 
keeping parameters. We will denote in particular the n/2-th 
PN order by 0(n), i.e. 0(n ) = 0(c~ n ). 


II. KHRONOMETRIC THEORIES WITH A MOND 
NON-RELATIVISTIC LIMIT 


The action of Hofava gravity ll29l[30 1 can be written as 


Sh — 


1-/9 

167tG 




+ 9u,v) > 


(1) 


where the spacetime has been foliated in spacelike hyper¬ 
surfaces, and the metric g )iv has been accordingly decom¬ 
posed in 3+1 form, i.e. we introduce the lapse function 
N = (—g 00 ) -1 / 2 , the shift 3-vector N i: = g oi , the induced 
3-metric 7 ^ (as well as its determinant 7 ) and the ex¬ 

trinsic curvature 

Kij = ^ {d tlij - DiNj - DjNi), (2) 

with D, denoting covariant derivatives relative to the geome¬ 
try of the spacelike hypersurfaces (i.e. /A, 77 ,. = 0). The mat¬ 
ter part of the action is instead represented by S m , where the 
matter fields ip couple to the covariant four-dimensional met¬ 
ric g liv , so as to enforce the weak equivalence principle and to 
confine Lorentz violations in the gravitational sector (at tree 
level) 1541 . The Lagrangian density L\ r j, is the most generic 
one at quadratic order in derivatives (up to total divergences), 
i.e. 

L kh = KijK” - + YZ~p* )R+ iZTp aial ’ (3) 


where a* = di In N, K = Y' 7 Kjj is the trace of the extrinsic 
curvature, and a, j3 and A are dimensionless free parameters. 
The parameter a regulates (among other things) the relation 
between the “bare” gravitational constant G appearing in the 
action and the “Newtonian” gravitational constant G y mea¬ 
sured by a Cavendish experiments, which turns out to be 


G n = 


2 G 

2 - a ‘ 


(4) 


The L 4 and L 6 Lagrangian densities are instead of fourth- 
and sixth-order, respectively, in the spatial derivatives /),, but 
contain no time derivatives H22ES3E3- This ensures that the 
theory does not suffer from the Ostrogradski instability E30 


3 See also section 2 of Ref. 03 for a pedagogical review of the Ostrogradski 
instability. 
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and most of all provides the anisotropic scaling necessary for 
power-counting renormalizability EH ED. For dimensional 
reasons, the L4 and L e terms must be suppressed by an en¬ 
ergy scale M*. This scale must be M* < 10 16 GeV to ensure 
that the theory remains perturbative at all scales, which is a 
necessary condition for power-counting renormalizability ar¬ 
guments to apply. Also, experimental constraints put lower 
bounds on M„. More precisely, to ensure agreement with 
submillimeter experiments ID El, it must be M* > 10 2 eV, 
and even more stringent bounds may be possible depending 
on the details of the percolation of the Lorentz violations in 
the matter sector beyond tree level. Indeed, observations of 
the synchrotron emission from the Crab Nebula show that this 
percolation should be suppressed if the theory is to remain vi¬ 
able and perturbative on all scales El- Several mechanisms 
have been proposed to suppress the percolation of Lorentz vi¬ 
olations from the gravity sector into the matter one, includ¬ 
ing fine tuning, “gravitational confinement” (58), “custodial 
symmetries” (e.g. softly broken supersymmetry ll59l l60l ), 
or dynamical emergence of Lorentz symmetry at low ener¬ 
gies in the matter sector, e.g. due to renormalization group 
flows [| 6 ll 1,62]. We refer the reader to Ref. ||63l for a re¬ 
view of these possibilities, and assume in this paper that one 
of these mechanisms suppresses the percolation to acceptable 
levels, so that the bound on M+ is M a b s < M* < 10 16 GeV, 
with M 0 b s ^ 10 ~ 2 eV. For these values of M* and at the 
low energies typically characterizing astrophysical observa¬ 
tions, the higher-order terms L 4 and Lq are typically negli¬ 
gible l49l . with the possible exception of black holes (whose 
causal structure does depend on the presence of the L 4 and 
Lq terms, c.f. the concept of universal horizon [[431 |46l ). 
When those terms are neglected, Horava gravity coincides 
with “khronometric” theory ll29l I30H . i.e. a theory with the 
action ( 0 , but with L4 and Lq set to zero ab initio. 

A useful way of writing the action of khronometric theory 
is to introduce a scalar field T (the “khronon”) defining the 
3+1 foliation, i.e. such that the constant-T surfaces coincide 
with the foliation’s spacelike hypersurfaces. Because of this 
requirement, this scalar field must have a timelike gradient, 
i.e. dj/fd,/!' < 0 within our conventions. In terms of this 
khronon field, the action of khronometric theory [i.e. Eq. 0 
with L 4 = Lq = 0] can be written in covariant form as |34„ 

ED 


ation, and 


p 

II 

<1 

3 

■g: 


(7) 

e = vx 

- 

(8) 

V {v n n) + a (/i n y) 

(9) 


(with 7 M1/ = g^ v + n^riv the projector onto the spacelike hy¬ 
persurfaces) are the acceleration, expansion and shear of the 
congruence defined by i.e. V M ri„ = —a„n p + < 7 M „ + 
§ 07 ^ 1 /. [Note that the vorticity + aun u ] = 

d[v n fi\ + a [n n v] vanishes identically because of Eq. 0.] 

It should be noted that the action 0 is very similar to that 
of Einstein-TEther theory IMl El, with the caveat that in that 
theory the vector is assumed to be timelike and unit-norm 
(thus = — 1 ) but not hypersurface orthogonal, i.e. n^ is 
a full-fledged (timelike and unit-norm) vector that cannot be 
expressed in terms of a scalar through Eq. 0 at the level of the 
action. For this reason, the vorticity of n M is not zero, and the 
most generic action for Linstein-/Lther theory is obtained by 
adding to the action 0 an extra term (c w being a 

dimensionless coupling constant), as well as a term + 

1) (where £ is a Lagrange multiplier) enforcing the unit-norm 
timelike character of the vector held n. 

Reference f50l proposed to modify the action of khrono¬ 
metric theory at the very large scales (i.e. very low energies) 
characterizing cosmological observations, i.e. in the infrared 
limit. The idea, as we outlined in the introduction, is that the 
cosmological evidence for Dark Matter comes from systems 
with accelerations a < ao ~ //<)/ 6 , and the theory intro¬ 
duced in Ref. Eq) seeks to reproduce the Dark-Matter phe¬ 
nomenology without any actual Dark Matter (with the pos¬ 
sible exception, as explained above, of some “dark baryons” 
on galaxy-cluster scales) by modifying the gravity theory in 
that low-acceleration regime. This corresponds to modifying 
the gravity theory on cosmological scales > 1 /ao, or equiva¬ 
lently energies < h ao ~ 10 34 eV. More precisely. Ref. Il50l 
considered a modified khronometric theory with action 


S = 


167tG . 


Jd 4 x^g R + f(a) + g^), (10) 


s =vtGh^[ R -l^ +3x)e2 


- panvcr^ 1 ' + aa^ 


T" 'S’mat( < Fj (5) 


where g is the metric’s determinant, R is the (four¬ 
dimensional) Ricci scalar. 


n u = - 




j-^dcFdpT 


( 6 ) 


is the (timelike) unit-norm vector field orthogonal to the foli- 


with a = y/j and n M still given by Eq. 0 . Refer¬ 
ence ll50l then showed that in order to obtain a MOND-like 
phenomenology in the non-relativistic, low-acceleration limit, 
the free function /(a) must asymptote to /(a) ~ — 2 Ao + 
2 a 2 — 4a 3 /(3ao) (where A 0 is a constant) for a<Cao, while 
they propose the limit /(a) ~ — 2A 0 b s (A 0 bs being the mea¬ 
sured cosmological constant) for a ao in order to reproduce 
GR (with a cosmological constant) in the high-acceleration 
regime. As we will show below, however, this theory does 
not produce a perturbative post-Newtonian (PN) expansion in 
time-dependent situations such as those of interest for cos¬ 
mology and astrophysics, i.e. the PN expansion turns out to 
be strongly coupled. We will show, however, that this prob¬ 
lem can be avoided with a slight modification of the theory of 
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Ref. f50ll . namely one with action 


S = 


107 tG . 


'dPxy/^g R— -(/? + 3A)6 > 2 - (3a 


f(a) + *5rnat(^: 9/ns)i ( 11 ) 


where again a = y/j^JaRa !7 , n fl is given by Eq. and /(a) 
satisfies again the asymptotic limit /(a) « — 2A 0 + 2a 2 — 
4a 3 /(3ao) for a <C ao. Note that this action can be rewritten 
in a 3+1 foliation adapted to the khronon, in the same way in 
which khronometric theory can be written in the two equiva¬ 
lent forms ([I}-(|3} and (|5j, thus obtaining 


s » = I - 1±4a- 

+ ^ _ p ^ g , ( 12 ) 

where a = ^ 7 ^ 0 +a" = ^/'yijtfa? in 3+1 form. In the high- 
acceleration regime relevant for astrophysical and experimen¬ 
tal tests (i.e. submillimeter, solar-system and pulsar ones), we 
impose that the theory reduces to khronometric gravity (plus 
a cosmological constant), i.e. for a ao (but a <C M*) we 
choose f(a) ~ —2A + aa 2 Q while at higher energies (i.e. 
a 2 > M*.) we may identify our theory with the full Horava 
theory. 

Of course, it remains to be seen whether the 
renormalization-group flow is compatible with this choice 
for the coupling function /(a), i.e. whether the MOND-like 
theory of Ref. fl50l (or a similar one, c.f. discussion in section 
|VI| ) is an infrared fixed point of the renormalization-group 
flow of Horava gravity. From this point of view, our treatment 
is purely phenomenological. 


III. THE HIGH-ACCELERATION REGIME 


As discussed above, for high accelerations (i.e. high en¬ 
ergies) a 7 > ao, our theory reduces to Horava gravity. In 
particular, for the accelerations ao < a < M* relevant for 
experiments on Earth and in the solar-system, as well as for 
most astrophysical (non-cosmological) observations, the the¬ 
ory described by actions ( [TT| or ( | 1 2 | reduces to khronometric 
theory. Here, we therefore review the experimental constraints 
on the coupling constants a, (3 and A of khronometric theory. 
Clearly, those constraints also apply to our theory. 

A linear expansion of the field equations of khronomet¬ 
ric theory on a Minkowski background shows that the theory 
presents a spin -2 graviton polarization propagating with speed 


4 A is related to the measured cosmological constant A,,| )s by A = 

Aob sG/Gc, where G c = 2G/(2 + ft + 3A) is the gravitational con¬ 

stant appearing in the Friedmann equations (64). In pratice, since BBN 

metal abundances and binary-pulsar observations constrain |/3, A| < a few 
x 0.01 14II1421 (c.f. also section |ul} , it must be A ~ A 0 b s - 
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P 

FIG. 1. Allowed parameter space for /? and A |4T1 [42l . from ob¬ 
servations with accelerations a > ao (solar-system tests, absence 
of Cerenkov vacuum radiation, BBN, pulsars) and from stability re¬ 
quirements (no gradient/ghost instabilities). 



Ct, as well as a spin-0 one with propagation speed c s . These 
speeds are given by EH go) 


r 2 = 

* 1 - p ’ 

2 _ (a — 2)(/3 + A) 

' s ~ a(/3 — 1)(2 + (3 + 3A) ' 


(13) 

(14) 


To avoid gradient instabilities on Minkowski space, one must 
impose c 2 > 0 and c? > 0. These conditions also ensure that 
energies are positive Eiisa, thus avoiding ghost instabili¬ 
ties. Even more stringently, to prevent ultra-high energy cos¬ 
mic rays from losing energy to gravitons by vacuum Cerenkov 
radiation 138) . the gravitational modes must also propagate lu- 
minally or superluminally, i.e. c 2 t> 1 and c 2 > 1 . 

To ensure that khronometric theory agrees with experi¬ 
ments at the level of the solar-system, one can solve the field 
equations at first PN order, and compute the Parametrized PN 
(PPN) parameters g). All these parameters turn out to be the 
same as in GR, with the exception of the preferred-frame pa¬ 
rameters aq and a 2 EMQ) 


cti 


CH 2 


4(a - 2/3) 

P-l 


(15) 


(cn — 2/3) [— (3(3 + /? + 3A) — A + ct( 1 + (3 + 2 A)] 


(a- 2)08- 1)00 +A) 


(16) 


Solar-system tests constrain |ai| < 10 4 * * and |a 2 | 
10” 7 0 . To satisfy these bounds, one can simply impose 
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a = 20 + 0(a\iOi 2 ) at leading order in a± and < 22 . This is 
sufficient to satisfy the constraints on both a.\ and a?, since 
both quantities are proportional to the combination a — 20. 
This allows decreasing the dimension of the theory’s parame¬ 
ter space from three (i.e. a, j3, A) to two (i.e. 0, A)rl 

Once the constraints discussed above are accounted for, the 
viable parameter space (0, A) is given by the cyan region in 
Fig. □ Additional bounds on the parameters then come, as 
mentioned above, from the requirement that BBN produce the 
observed element abundances f4Tll42l l64l l 66 l (orange region 
in Fig. [TJ. Also, stringent bounds (represented in green in 
Fig. 0 come from the absence of any anomalous precession 
in observations of isolated pulsars mm, as well as from 
the change of the measured period of binary pulsars under 
gravitational-wave emission E] [42]. Indeed, the latter ef¬ 
fect puts very strong constraints on (3 and A, because the pres¬ 
ence of a khronon field coupled non-minimally to the metric 
causes the appearance of dipolar fluxes in the gravitational- 
wave emission from binary systems, besides the quadrupolar 
fluxes of GR ED 02- Because binary-pulsar observations 
are in good agreement with the GR predictions, these dipolar 
fluxes must be suppressed by sufficiently small values of the 
coupling constants. 

Nevertheless, as is clear from Fig. [TJ there is a sizeable re¬ 
gion of parameter space where khronometric theory [and thus 
the theory described by Eqs. © or is viable, around 
the limit (3 = A = 0 (in which GR is recovered at high 
accelerations). Note that in this viable region of parameter 
space, black-hole solutions that arise from gravitational col¬ 
lapse m have also been shown to exist EM!). These 
solutions present properties compatible with current electro¬ 
magnetic observations of black-hole candidates (i.e. their ex¬ 
terior geometry is very close to the black-hole solutions of 
GR) E2|49). 

Finally, as discussed in the previous section, our theory re¬ 
duces to Horava gravity in the UV regime a M*. There¬ 
fore, constraints coming from sub-millimeter tests of the 1 /r 2 
decay of the Newtonian attraction force are satisfied provided 
that M* > 10 ~ 2 eV 11321 , while tests of Lorentz invariance in 
the matter sector will be passed provided that a suitable mech¬ 
anism exists that suppresses the percolation of Lorentz viola¬ 
tions from gravity to the matter sector (c.f. discussion and 
references above). In addition, as alluded above, the higher- 
order derivative terms of Horava gravity are important for the 
propagation of signals in a black-hole spacetime, but are not 
expected to destroy its causal structure (which still possesses a 
universal horizon from which no signals can escape, not even 
with infinite propagation speed B31I461 ). 


5 One may also impose the bounds a 1 1 < 10“ 4 and 0:2 < 10 -7 without 
exploiting the fact that both a\ and a 2 are proportional to a — 2/3. How¬ 
ever, this would produce a one-dimensional parameter space, which turns 
out to be a subset of the two-dimensional parameter space that one obtains 
by choosing a = 2(3 + 0(a 1 , 0 : 2 ). In this sense one may choose a & 2/3 
without loss of generality, c.f. Ref. ED for a detailed discussion. 


IV. THE LOW-ACCELERATION REGIME 

In this section, we will study the 1PN expansion of the the¬ 
ory described by actions © or ©• While our treatment is 
valid in both the high- and low-acceleration regimes, we will 
focus mostly on the latter. Indeed, as discussed in the previous 
section, at high accelerations the theory reduces to khrono¬ 
metric theory /Horava gravity, for which the 1PN expansion 
has already been derived in Refs. If33l[40 l. and shown to agree 
with experimental constraints coming from solar-system tests 
in large portions of the parameter space. As a check of our cal¬ 
culation, we will however verify that we reproduce the 1PN 
expansion of khronometric theory derived in Refs. 63 ED, 
confirming in particular their expressions for the preferred- 
frame parameters a\ and <22 [Eqs. © and ©]. 

The calculation of the 1PN expansion in the low- 
acceleration regime, which we present below, may a priori 
be expected to be of purely academic interest. After all, the 
tests of the PN dynamics of GR (solar-system tests and binary 
pulsars) are in high-acceleration regimes, while systems with 
accelerations a <C do (such as those encountered in cosmol¬ 
ogy) have velocities too small relative to the speed of light to 
test even the first PN order (with the accuracy of current data). 

Nevertheless, we will show that surprises arise in the course 
of the calculation. In particular, we will show that if one sets 
the couplings (3 and A to zero (as in the original theory of 
Ref. f50l ) or to values below a certain threshold, the 1PN 
expansion in the low-acceleration regime becomes strongly 
coupled. We will show that this prevents the theory from re¬ 
producing the Dark-Matter phenomenology at accelerations 
a -C do, at least in dynamical/time-dependent situations and 
unless the couplings j3 and A are significantly different from 
zero. We will show this explicitly by calculating the rotation 
curves of the gas surrounding a galaxy (whose mass grows due 
to accretion from the intergalactic medium - IGM) at 1PN or¬ 
der. Based on this toy model, we will then compute a lower 
bound on the combination |A + 0\ [c.f. Eq. f73j)], i.e. we will 
determine the minimum value of this combination for which 
the theory avoids the aforementioned strong-coupling prob¬ 
lem, and can thus reproduce the Dark-Matter phenomenol¬ 
ogy at low accelerations. We will show that by combining 
this bound with existing constraints on the couplings from the 
high-acceleration regime the theory remains viable in a non- 
negligible region of parameter space (c.f. Fig. [2]). 


A. Modified field equations 


As a first step toward computing the 1PN expansion, let us 
first derive the field equations by varying the action in adapted 
coordinates, i.e. Eq. ©. The variation with respect to the 
lapse N gives 


( 3 ) i? 


K 1 ' 3 K l? + 


1 + A 

1-/3 


K 2 + 


/(«) 

1-/3 


2 

1^0 


\a 


2 

1-/3 


A(xa‘) 


16irG£ 


(1 - /3)c 4 ’ 

(17) 
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the variation with respect to the shift iV, gives 


1), K IJ - 


1 + A 
1-/3 


7 ij K = - 


8 trGJ l 
(1 ~P)c 


(18) 


and the variation with respect to the 3-metric 77 yields 


1 


1-/3 


{,) R ij (k ij k 


+ ^ D k _ A -l±| 7 i)fc)^ + 2K ik K° k 


1 + /3 + 2A 
1-/3 
1 

(l-/3)/V 
1 


1 


K VJ K - 77 ^ ( + 


(D l D J N - 7 ij D k D k N) 


1 + A 
1-/3 


JG 


1-/3 


1 7 

X« a 


/( a ) 

2(1-/3) 


7*J = 


8 nG 


(l-P)c 


rV 3 . 


(19) 


In these equations, \ = f (a)/(2a), Di denotes the covariant 
derivative compatible with 77 , while is a shortcut for d t — 
N k D k . Also, the terms £. J'\ T' :> come from the variation of 
the matter action, i.e. 


similar to GR, i.e. Eq. ( fl9| is a modified evolution equation 
and Eq. ( fT 8 f > is the modified momentum constraint. On the 
other hand, Eq. ( fi~7| ) clearly looks like a modified Hamiltonian 
constraint, but a key difference from GR is present. Indeed, in 
GR one may in principle choose a specific gauge (defined by 
some conditions on N and Ni), choose initial data compatible 
with the constraints, and evolve the evolution equation, which 
would ensure that the constraints are satisfied at later times. 
This is not possible in the case of Eqs. ([T7|)-([T9|), since we 
have already used up our “time” gauge freedom by adapting 
our coordinates to the preferred foliation. This can be seen 
explicitly by transforming the action of Eq. © to that of 
Eq. ( fl2| ), by choosing a 3+1 decomposition such that the time 
coordinate t matches the khronon scalar T. As result, the lapse 
N is not a gauge field in Eqs. ([T7|)-([T9|, but should rather be 
solved for at each step of the evolution via Eq. ©• Indeed, 
it can be shown that once Eqs. ( p~ 8 j ), ( [T9 | i and the equations of 
motion of matter are assumed to hold, Eq. ( fT7) ) is needed to 
ensure the validity of the khronon evolution equation (which 
is obtained by varying the covariant action ( [TT| with respect 
to T) USED. Also, as we will see below, the lack of free¬ 
dom to “gauge away” the lapse will be the origin of the PN 
strong-coupling problem mentioned aboveQ 


£ = - 

1 f^mat 

( 20 ) 

\fl SN ’ 

1 ss mat 

T = 

( 21 ) 

Vi sn z ’ 

'J'ij — 

2 <5S mat 

( 22 ) 

N Vi hij 


and are related to the canonical stress-energy tensor compo¬ 
nents, = ( 2 /yf^g) {SS^/Sg^), by 


B. Post-Newtonian expansion 

To calculate the PN, let us start by writing the most generic 
perturbed flat metric in Cartesian coordinates (x° = ct,x z ) 
(see e.g. Refs. HI1): 

3oo = -1 - 4</> - \<t> m + 0(6) 


£ = N 2 T 00 , 

(23) 

90i — - 3 H- 

C J c, 

/ 2 ' 

Ji = N ( T 0i + N i T 00 ) ; 

(24) 

9a = y 1 - ^ 

-J~ij _ jiij _ 00 

(25) 

i ^ > 


+ + 


Cij 


0(4). 


(27) 


(obtained by contracting that equation with 77 ), we obtain 


S) R 2 / 1 + A\ ,, 

- — ( 1 - 3--- ) D t K + 3K u K kl 


1-/3 N \ 1-/3 

1 + 2/3 + 3A 2 4 

H-:---/V — 


1-/3 


1-/3 


Xa 2 = 


N{ 1-/3) 

167tG 


D k D k N + 


1-/3 


/(«) 


(1 — /3)c 


T+-N k J k ), (26) 


Under transformations of the spatial coordinates, 
w, 0 , </> (2) transform as scalars, vj t . (i behave instead 
as transverse vectors (i.e. diW 1 = diQ = 0 ), and £7 is a 
transverse and traceless tensor (i.e. d^ = 0 ). 

Since we have already chosen our time coordinate to coin¬ 
cide with the khronon field T, we only have freedom to rede¬ 
fine the spatial coordinates on our foliation, i.e. we are only al¬ 
lowed to perform gauge transformations -+ h /JiV + d( ll ^, 


which will come in handy later. 

Several comments are in order about these field equations. 
First, for A = /3 = 0 they reduce to those presented in 
Ref. ff5Qj[ f] Also, the structure of these equations is clearly 


6 Note that our definition of f(a) as given in the action [Eqs. ED or 0i 
differs by a factor —2 from the definition chosen in Ref. ED. 


7 Although in the next sections we will show that the PN dynamics becomes 
non-perturbative (i.e. strongly coupled) when j3 and A are equal or close 
to zero, this does not necessarily mean that the theory’s structure itself is 
pathological, even when (3 = X = 0. Indeed, one may in principle inte¬ 
grate in time the evolution equation 0 as in GR, and solve the Hamilto¬ 
nian constraint 0 (which is an elliptic equation for N) at each time-step 
(given appropriate boundary conditions). This would be possible even for 
/3 = A = 0, although the resulting dynamics would not be perturbatively 
close to the Newtonian one. 
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with representing the perturbation and = 

(0 ,&) a purely spatial vector. For this calculation, we find it 
convenient to impose the gauge conditions £ = Q = 0 f 68 l . 
As a result, the lapse, shift, spatial metric and acceleration at 
1PN order are given by 


mass distribution of matter is different than in GR). This is 
important as it allows the theory to reproduce the successes of 
the ACDM model in the interpretation of gravitational lens- 
ing from galaxies and clusters of galaxies |[T9ll50ll70S . Also, 
based on Eq. ([33]), we can write 


N = 




,oo 


— 1 + — 


AT 1 i 

Ni — got — T 


7 ij =9ij= 1 - 


di4> <pdii 

a i = — ~ 2 — 

C A c 4 



(28) 

■£+<*»). 

(29) 

) ’ 

(30) 

+ dl ^ + 0 ( 6 ) 

(31) 


i/j = cl ) + S 4-+0( 4), (34) 

C“ 

where we have defined the potential Sip, which will appear in 
the rest of the calculation [c.f. Eq. ( |4 1 [ i below]. 

Using this result in Eq. © to lowest order in 1 /c we ob¬ 
tain l50l 



AttG P + 0(2), 


(35) 


which can be used to compute the left-hand sides of the field 
equations To compute the right-hand side of those 

equations, we use a perfect fluid stress-energy tensor, i.e. 

T^=(p+^ju^+pg^, (32) 


where since /(a) oc a 2 = 0(4), one has that x = f'(o) /(2a) 
is of zeroth-order in 1/c. In the high-acceleration regime, 
/(a) ss «a 2 , thus \ = a and this equation becomes the usual 
Poisson equation 

V 2 <p N =4nG N p + 0{2), (36) 


where p is the matter mass-energy density, p the pressure and 
wT = dx 1 */dr the four-velocity of the fluid elements (with r 
the proper time). 

Before proceeding with the calculation, let us clarify the 
PN order of the function /(a). As mentioned in section [5] 
in the high-acceleration regime (i.e. for ac 2 ao), /(a) ~ 
—2A + aa 2 ss aa 2 , so Eq. © implies /(a) = 0(4). Note 
that in deriving this scaling we have used the fact that A is 
comparable to the observed value of the cosmological con¬ 
stant, i.e. c 4 A ~ c 2 Hq ~ a§ -C a 2 c 4 , which allows neglect¬ 
ing the —2A term. (Of course, this corresponds to the known 
fact that the cosmological constant has negligible impact on 
the 1PN dynamics on small scales.) In the low-acceleration 
regime (i.e. on cosmological scales), the cosmological con¬ 
stant would instead be expected to enter the 1PN dynamics. 
Indeed, for ac 2 <C ao, /(a) « — 2A 0 + 2a 2 — 4a 3 c 2 /(3ao), 
and assuming (as is natural to do) that Ao is comparable to the 
observed value of the cosmological constant, the term — 2 Ao 
dominates over 2a 2 — 4a 3 c 2 /(3ao) = 0(4). However, in 
order to have the same scaling /(a) = 0(4) as in the high- 
acceleration regime, we can simply move the cosmological 
constant to the right-hand side of the field equations, and ab¬ 
sorb it in the matter stress-energy tensor as a “fluid” com¬ 
ponent with equation of state p/c 2 = —p = —Ac 2 /(8nG), 
as routinely done in cosmology. Therefore, in what fol¬ 
lows we will consider /(a) = 0(4) in both the high- and 
low-acceleration regimes, with the caveat that in the latter 
/(a) w 2a 2 — 4a 3 c 2 /(3ao) and the matter is meant to include 
a “Dark-Energy” component p/c 2 = — p = —Ac 2 /( 8 ttG). 

With these Ansatze and scalings, deriving the 1PN field 
equations is now straightforward. In particular, expanding 
Eq. ( |26) to lowest order in 1/c yields ll50l 

ip = <t> + 0(2), (33) 

which implies that light deflection behaves as in GR (except, 
as we will show below, that the relation between <p and the 


with Gn given by Eq. 0 . At intermediate and low accelera¬ 
tions, x is not necessarily constant, and defining an “interpo¬ 
lation function” 


fi = 1 - |, (37) 

Eq. ([35) becomes the modified Poisson equation of the 
MOND dynamics fl9ll22ti24l. i.e. 



4nGp + 0(2). 


(38) 


In particular, in the low-acceleration regime ac 2 <C ao (i.e. in 
the “deep-MOND regime”), /(a) « 2a 2 — 4a 3 c 2 /(3ao) and 
Eq. ( |35) becomes 

^V<A=47rGp. (39) 

ao 



From the off-diagonal part of the modified evolution equa¬ 
tion [Eq. ( [19) ] we obtain 

Cij = 0 ( 2 ), (40) 

i.e. appears at higher order than 1PN. This is of course ex¬ 
pected, since this term represents gravitational waves, which 
do not enter in the 1PN metric in GR. 

Solving then the trace of the evolution equation [Eq. ©] 
to 0(4), we obtain 

-^V 2 Ai p = + -^(-24irGp-8npv 2 - TdiCpdiCp 

c 1 2 \ 

+ x^itp^itp — 8 (pV~<p + (2 + P + 3A)(<9 t V 2 w + 8 d 2 (j .>)^ , 

( 41 ) 
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and by replacing this expression in the modified Hamiltonian 
constraint [Eq. ( fl7| )] at 0(4) we find 


V • 



47rGp + c 2 


/O) 

2 


+ 4 f 8 tt Gpv 2 + 12t rGp + 2V0 • V</> - ^xW • V</> 
c z V 2 

—-(2 + ^ + 3A)(9tV 2 w + 39 2 ^)^. (42) 


Finally, the 1PN equation for the “frame-dragging” potential 
Wi can be obtained from the momentum constraint [Eq. (p~8|)], 
whose expansion yields 

v 2 «ii + 2 (|il)a.v 2 ^ 

= V»Gpv, _ 2 ( 2+/ + ^ ) W . 

1-/3 V /3 — 1 / ,tv 

(43) 


By taking the divergence of this equation we obtain 

V 2 'V 2 w = —J— [87rG9 t p - (2 + /3 + 3A)5 t V 2 ^], (44) 
p + A 


where we have used the condition diW 1 = 0 (c.f. the defini¬ 
tion of Wi), and the energy conservation to Newtonian order, 
d t p = —d i (pv' l )[ 1 + 0(2)]. Denoting by <j> N = 47 rGjvV _2 p 
[with Gn given by Eq. (|4])] the Newtonian potential in the 
high-acceleration regime, we can rewrite Eq. ( |44| in the more 
useful form 

V 2 w= -^d t [(2-a)<^v-(2 + /3 + 3A)0]. (45) 

Note that Eqs. (|42|-(|45]i are valid both in the high- 
acceleration regime, in which /(a) ss aa 2 , and in the low- 
acceleration, deep-MOND regime, characterized by /(a) ss 
2a 2 — 4a 3 c 2 /(3ao) (and thus x = 2 — 2c 2 a/ao). In the high- 
acceleration regime, we must of course recover the known 
results for khronometric theory, namely that all the PPN pa¬ 
rameters vanish except for ax and a 2 , which are given by 
Eqs. ( p~5j ) and > [T 6 | ). Indeed, we show explicitly that this is 
the case in the Appendix. 

The low-acceleration, deep-MOND regime is instead ana¬ 
lyzed in detail in the next section. However, already looking 
at Eq. m we can understand that the 1PN expansion in the 
deep-MOND regime may have a non-perturbative character, 
because the right-hand side seems to diverge for /3 + A — > 0. 
Clearly, this cannot be the case in GR, where we know that 
the 1PN expansion is perturbative. Indeed, in GR one has 
a = /3 = A = 0 and (j> = cf>N, thus the two terms in round 
brackets on the right-hand side cancel out. This is consis¬ 
tent with the fact that in GR one can set w = 0 by a gauge 
transformation of the time coordinate |[ 68 l (while still impos¬ 
ing the conditions £ = Q = 0 by a gauge transformation of 
the spatial coordinates, as we do in this paper). Because in 
the khronometric theories that we are considering we already 
fixed the time foliation by adapting it to the khronon T, we 


have no residual gauge freedom to set lo to zero, and V 2 w 
may indeed diverge (in general) when A, /3 —» 0. Another 
way of seeing that the case /3 = A = 0 is pathological is to 
note that if we had started from such a theory, we would have 
derived Eq. ( |43| with /3 = A = 0, i.e. the same equation as in 
GR. That equation, however, would have no dependence on lo, 
which would therefore remain completely undetermined. This 
is not a problem in GR, as w is a gauge mode (so it should in¬ 
deed remain undetermined), but is a problem in the modified 
khronometric theory of Ref. Il50i . because lo is not a gauge 
mode there. 

Indeed, already in the high-acceleration regime the terms 
in round brackets on the right-hand side of Eq. ( |45j ) do not 
cancel out (in general) in the theories we are considering. 
This is because (f> = is that regime, but a, /3 and A 
are in general non-zero. Of course this corresponds to the 
fact that in khronometric theory the preferred-frame parame¬ 
ter 0 L 2 becomes large when A + /3 is small, unless a ss 2/3 
[c.f. Eq. ©]0 For high accelerations, however, we have 
already discussed that one does indeed have the freedom to 
set a ~ 2/3, so as to satisfy the solar-system constraints 
|ai| < 10 - 4 and |qt 2 | 10 '. There is therefore no strong¬ 

coupling problem in the viable part of the parameter space of 
the couplings at high accelerations. 

The situation is different in the low-acceleration, deep- 
MOND regime, since <p ^ <pN there. Indeed, we will show ex¬ 
plicitly that the right-hand side of Eq. ( |45] > diverges in the limit 
/3, A —> 0, in low-acceleration, time-dependent/dynamical 
systems. We will also show that this strong-coupling problem 
appears in a region of parameter space that would be other¬ 
wise allowed based on experimental constraints coming from 
the high-acceleration regime. 


C. The strong-coupling problem in galactic rotation curves 


As shown in the end of the previous section, the 1PN equa¬ 
tions present a strong-coupling problem in time-dependent sit¬ 
uations, if A and /3 are very close to 0. In this section we will 
show this explicitly by solving the 1PN equations for a toy 
model consisting of a spherical galaxy whose mass M(t) in¬ 
creases linearly as a function of time due to e.g. accretion of 
gas from the IGM. We will then compute the conditions that A 
and /3 must satisfy to avoid this strong-coupling problem, and 
show the resulting parameter space in which the theory re¬ 
mains viable. More specifically, we will compute the rotation 
curves for such an accreting galaxy outside its luminous cen¬ 
ter, assuming that no Dark Matter is present, and assess for 
what values of A and /3 the aforementioned strong-coupling 
problem modifies the rotation curves in a way that is incom¬ 
patible with observations EMI. 


Indeed, for a = 2(3 and (j> = </>n, the right-hand side of Eq. (45} is 
independent of /3, A. 
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1. The Newtonian order 


At Newtonian order, the equation for the perturbation <fi 
at low accelerations |V^>| -C a o is given by Eq. ( |39| ). In 
spherical symmetry, however, (j) is only a function of the dis¬ 
tance r from the galaxy’s center, thus we can always write 
|V0|V<^>/ao = VS 1 for some scalar function S(r). Insert¬ 
ing this definition in Eq. |39|, we obtain that S must coincide 
with the GR Newtonian potential cj)^. Therefore, to find the 
MOND gravitational potential <fi in spherical symmetry, we 
can simply solve the corresponding Newtonian problem in GR 
for and then compute (!) by solving 


d <f>(r) 
dr 



#jy(r) 

dr 


(46) 


As our toy model for an accreting galaxy, let us consider a 
spherical body with mass M = Mq + Mt (with M and Mo 
constants) and radius R, surrounded by a spherically symmet¬ 
ric, stationary accretion flow [whose density, simply by mass 
conservation, is p = M/[4irr 2 v r (r)\, where v r (r) is the ra¬ 
dial infall velocity as a function of radius]. Note that because 
no Dark Matter is assumed to exist, we identify R with the 
galaxy’s half-mass radius, which is related to the (baryonic) 
mass M by the observational fit fZ 3 


'-5.54 +0.56 log 10 (^) 
for iogio (^) > 10.3, 


logi 0 (i?eff/kpc) = < 


-1.21 + 0.14 log 10 (^) 
for iogio (]|t) - 9 10 ' 3 - 

(47) 

Let us focus on the region outside the galaxy’s radius R, 
where only the accreting gas and the cosmological constant 
are present. In this region, cf>N is given by 


<t>N =-— + 0finite(M, A obs ) , (48) 

r 


At distances from the body’s center r«r o (but r > R, i.e. 
outside the galaxy), (j) coincides with as given by Eq. ( |48| , 
while for r !$> r 0 , <f> is given by Eq. <(46]». We can therefore 
assume a sharp transition at r = ro, and solve Eq. ( [46] ) by 
imposing continuity of and its first derivative, i.e. 4 >(ro) = 
~G N M/r 0 + G>fi ni te(M, A obs ) and d(j)/dr = G N M/rl + 
C’finite(M, A obs ), thus obtaining 

cj) = \J GVMao ^ln ^ — 1^ + C7flnite(Af, A obs ) (50) 


for r > r 0 . 

Finally, as we will show explicitly in the next section, we 
do not need the explicit form of the potential 4 > = (j)^ inside 
the galaxy (i.e. for r < R <C ro) to solve the 1PN equations, 
if we focus on the terms that dominate when /3 + A —> 0. 


2. The metric at 1PN order 


At 1PN order, the metric is characterized by the potentials 
w, w l and <j> (2) . In spherical symmetry, however, w % = oQ 
To determine w, let us start from Eq. By using the 

Green function of the Laplace operator we obtain 


u)(x, t ) = 

-*[/, 

\_J r'>r o 


d 3 £ ( 2 - - (2 + /3 + 3A)<j5(r', t) 


47 t(/3 + A)|x — x'\ 


f d^x' @ 3X)(t>N(r\ t) 

I r' <ro X 4?r(^+ A)|f - X'\ 


+ ^o, (51) 


where +o is an integration constant, and we have used the fact 
that 4> = (f’N at high accelerations (i.e. for r < ro). As 
already noted in the previous section, if a ~ 2/3 (as required 
by solar-system tests), the second integral in Eq. is finite 
when /3, A —>• 0, and can therefore be neglected with respect 
to the first one, which diverges. More precisely, by assuming 
a = 2/3 + 0 (ai,a 2 ) (so as to pass solar-system tests), in 
spherical coordinates the above solution becomes 


where Ofi n i te (AI, A obs ) denotes corrections (proportional to 
either A obs or M) that remain finite as /3, A —► 0. Indeed, these 
corrections are clearly independent of /3, A in this case, and 
are also time-independent, because both p = A obs c 2 /( 87 tGat) 
and p = M/[Anr 2 v r (r)\ do not change with time. 

To compute <j>, one can then just solve Eq. ( |46[ >. To do so, 
one needs to specify conditions ensuring a smooth transition 
to the GR solution, which is valid in the high-acceleration 
regime near the galaxy. In particular, let us define the tran¬ 
sition radius 


t o 


G]yM 

a 0 


(49) 


w(r,f) 


1 

W+ A) X 


d t {- I drV 2 [2(l - P)ct) N {r',t) - (2 +/3 + 3A)c/>(r',f)] 
l r J r 0 

+ J drV[ 2 (l - f})(j) N {r',t) — (2 +/3 + 3A)^(r', f)]| x 
x [1 + 0(a i, 0 : 2 )] + f/’o + G(/3 + A ) 0 . (52) 


Because 0 diverges as In r as r —> 00 , the second integral on 
the right-hand side of this equation formally diverges. This 


at which the Newtonian gravitational acceleration |V0A/j = 
GM/r 2 + 0 finite (AT, A obs ) matches the acceleration constant 
do- [Note that ro is larger than the half-light radius given by 
Eq. (|47j) for typical galaxy masses.] 


9 This follows from the requirement that diW 1 = 0. imposing regularity at 
r = 0. Alternatively, one can solve the divergenceless part of Eq. ID- 

noting that in spherical symmetry the velocity only has a radial component 
'iv(r), which can be expressed as the gradient of a scalar potential. 
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is simply because the PN formalism is by definition a pertur¬ 
bative expansion on a Minkowski background [c.f. Eq. m- 
Of course, for any given spacetime one can choose locally 
Riemannian coordinates x a centered on a given event, and 
such that the metric is locally = 77 y v + 0(r/lZ) 2 , where 
r ss \Jr\ a px a xP is the proper distance from the event and 
7Z is the curvature radius of the spacetime at the event. In the 
particular case of a system embedded in a cosmological space- 
time, 7Z ~ c/ H ( H being the Hubble rate), i.e. the Minkowski 
metric is the appropriate background metric only on length- 
and time-scales much smaller than the cosmological ones (i.e., 
respectively, the Hubble radius and Hubble time) 01. For this 
reason, we can truncate the second integral on the right-hand 
side of Eq. ( [52] ) at a cut-off radius A max , which can be thought 
of as much smaller than the present Hubble radius but much 
larger than the typical size of the luminous component of a 
galaxy. 

In practice, the cut-off /? niax never enters our calculations 
and results, as it can be renormalized in the integration con¬ 
stant i/jq. Indeed, once this cut-off is imposed, we can use 
Eqs. ( [48] ) and ( |50| > for 0n and (f> in Eq. ( |52| ), and the integra¬ 
tion yields the following expression 


w(r,f) = 


M 


I ao Gn 


72r(/3 + A) 
— (2 + f3 + 3A) 
MRmax 


72G N (1~P) (r 2 + r 2 ) 


17r 3 + 28rg - 6 r 3 In ( — 


8(/3 + A) 


, M{t) 

— 16Gat(1 — p ) + i? max (2 + p + 3A) x 


cloGn 

M 


(5 + 21n(r 0 /A max ) 


+ Otinite(M, Acbs) + 0(13 + A) c 


x [1 + G(ai,a 2 )\ + i>o 
(53) 


from which it is clear that the terms that depend on the cut-off 
radius can be absorbed in the integration constant ipo. There¬ 
fore, the final solution for the potential w(r, t) is simply 


= -- 


M 


72r(/3 + A) 
— (2 + p + 3A) 


72G JV (l-/3)(r 2 + rg) 


ooGat 

~W) 


17r 3 + 28r 3 - 6 r 3 In ( — 
ro 


x [1 + 0(a 1 , 0 : 2 )] + Gfinit e(M, A 0 b s ) + G(/3 + A) 0 , (54) 


Let us now consider the modified Hamiltonian constraint 
Eq. ©. Because of spherical symmetry and taking into ac¬ 
count only the terms that diverge when (3 + A —► 0, at 1PN 


order that equation becomes 

2 d f. 2 d0 d<l> m \ 

a^r 2 dr V dr dr ] 

= _ 2 + g + 3A a, v ^ + 0()3 + A) ° = 

O 1 ft 1 Q \ 

- 2 ^ + A) dt [2(1 - P)<!>n ~ (2 + P + 3A)0] x 
x [1 + 0(a.i, 012 )] + G(/3 + A) 0 (55) 


By inserting the explicit expression for <j> at low accelerations 
[Eq. ( [50] )] and isolating the derivatives of (f> m on the left-hand 
side, in the deep-MOND region r > r 0 , Eq. (|55j) becomes 


d_ ( dPa, \ 
dr y dr J 


F(r,t) = 


~ { 2 t(pt\)ro 2 ^ [2(1 “ - ( 2 + P + 3A )^ X 

x [1 + 0(a 1 , a 2 )] + 0(p + A) 0 + Gfi n ite(A7, A 0 bs) (56) 


where F(r,t) represents the source on the right-hand side. 
To solve this equation, let us construct the Green function 
G(r, r'), i.e. the solution to d r (rd r G) = S(r — r'). As 
usual, the Green function can be constructed from solutions 
of the homogeneous problem. In brief, for r ? r', the equa¬ 
tion defining the Green function becomes d r (rd r G) = 0, 
which has the general solution G(r,r') = K 2 ln(r/ro) + 
Ki, with I\i and AT 2 being integration constants. Imposing 
then the junction conditions G| r=r / .. 0 - = G| r=r /_ 0 + and 
rd r G\ r - r : +0 + — rd r G\ r — r i +0 - = 1 to account for the pres¬ 
ence of the Dirac delta on the right-hand side, we then obtain 


G(r, r') = 


I< 2 In y—J + 
ln (?) + K 2 ln 


r < r' 


+ Ki r > r' 


(57) 


The general solution to Eq. ( |56| l can then be written as 

P Rmdx 


<t>m{r,t) = 


K ‘2 In ( - ) + A'i 

ro 


dr 1 ln ^ —^ F {?, t) 


dr'F {?, t) 


(58) 


which explicitly gives 

. . _ «o(2 + P + 3A) 2 Af 2 

<P(2>lr ’ ’ ~ 432 (p + A )M 2 

+ 12Jvi(rQ - A^ ax ) + 9Aii? 3 ax ln 


5(r 3 ^r 3 ) 


r 0 


3 ln 

\r 0 

+ 3I< 2 Ftf ln 


r 3 + A? + 4A' 2 (rQ - A 3 ax ) 
R„ 


u max 

ro 


x [1 + 0(a 1 , a 2 )] 


+ 0(P + A) 0 + Gfi nite (A7, A obs ). 


(59) 
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The integration constants K\ and K 2 can be fixed by im¬ 
posing that the solution matches the high-acceleration solu¬ 
tion </>pJ sh acc and its derivative at r = r$, i.e. 0 (2j (ro) = 

0ra gh acc (r 0 ) = Hi and d r (j> m (r 0 ) = d r 4>m Sh acc (r 0 ) = H 2 , 

thus obtaining 


4> (2) {r,t) = Hi + H 2 r 0 hi 



a o (2 + 0 + 3A) 2 M 2 
432(0 + A )M 2 X 


x ^ 5 ( r o ~ r3 ) + 3 (r 3 + 4r o) ln j j x 
x [1 + 0(a i, 02 )] + ©(0 + A) 0 + ©finite (M, A obs ). (60) 


Note that the cut-off radius f? max is once again absorbed in 
the integrations constants Hi and H 2 (as for the potential u> 
earlier in this section). Since in the high-acceleration regime 
the theories that we consider reduce to khronometric theory, in 
which no strong-coupling problem is present when 0 + A —>■ 0 
if the s olar-sy stem constraints are satisfied (c.f. discussion in 
section IV B 1 , we have Hi = ©(0+A)° and H 2 = ©(0-fA)°, 
and Eq. (|60|) can be rewritten simply as 

a ( rt )- a ° ( 2 + Z 3 + 3A ) 2A ’ /2 x 

’ ’ 432(0 + A)M 2 

x ^ 5 ( r o “ f 3 ) + 3(r 3 + 4 t"o) ln 

x [1 + 0(a i, a 2 )] + 0(/3 + A)° + ©fi n ite(Af, A 0 b s ). (61) 


3. The impact of the strong coupling on the rotation curves 


Let us now explore the impact of the strong-coupling prob¬ 
lem described above on the rotation curves of galaxies. To this 
purpose, let us model the gas (whose velocity is measured to 
determine the rotation curves) by test particles following cir¬ 
cular geodesics in the deep-MOND region r > r 0 . Because of 
spherical symmetry, we can assume that the orbits are on the 
equatorial plane, without loss of generality, i.e. in spherical 
coordinates x M = (cf, r, 0, tp) the four-velocity of the gas is 


u>i = ~T (c, °, °, 0), (62) 

dr 

where at 1PN order the relation between coordinate time t and 
proper time r is given by 


dt 

dr 


= 1 - 


</> 


0<i) 2 

2c 2 


+ 0(4), 


(63) 


which follows from the normalization condition u ll u 1 ' = —c 2 . 
The geodesics equation 


Focusing on the radial component, and because d 2 r/dt 2 = 
dr/dt = 0 for circular orbits, Eq. ([65]) then gives 


^(^,1PN 


2-2 90 r 

= r tp = r — [t,r ) + — 
or c z 


+ 2 0(f , r)~(t,r) + ^ (f , r) 


<90 

dr 
8 2 uj 
dtdr 


(t,r) 

(t,r) 


■0(4). 

( 66 ) 


By using Eq. ( |50| ), at the lowest (i.e. Newtonian) order, this 
equation yields 


< N = VG N Ma 0 + ©finite(M, A obs ) + 0(2), (67) 

i.e. the rotation curves of galaxies are flat in the deep-MOND 
region. [Note also that the scaling of Eq. © with the mass 
agrees with the Tully-Fisher relation for disk galaxies and the 
Faber-Jackson relation for elliptical galaxies and clusters, c.f. 
Ref. |[T9l for a review of these two relations in the context of 
MOND.] At 1PN order, and focusing on the terms that diverge 
as 0 + A —» 0, the rotational velocity becomes 

= VG N Ma 0 + ^(^^{t,r) + 

+ ©(0 + A) 0 + ©finite (A^) Aobs) + 0(4) (68) 


or more explicitly, by using the solutions given by Eqs. 
and ( [60] ), 


v:„ , m = \/G N M(t)a 0 


y 2 

-'<£,1PN 


1 r ao(2 + 0 + 3A) 2 M 2 


144(0 + A) M(t) 2 _ 


4 ( r o ~ r3 ) + 3r 3 ln — 


r o 


M 2 


I UqGN 


36r(0 + A)Af(f) Y M(t) 


(2 + 0 + 3A) x 


x ^4r 3 + 14rg — 3r 3 ln ^— J J + 36(0 — l)rj 

x [1 + ©(ai, at 2 )] + ©(0 + A) 0 + ©fi n ite(Af, A 0 bs) + ©( 4 ) • 

(69) 


Clearly, if M 0 and 0 + A —» 0, the 1PN terms in this 
expression will dominate over the Newtonian ones, spoiling 
the agreement with galaxy rotation curves and with the Tully- 
Fisher and Faber-Jackson relations. In the next section we will 
determine exactly for what values of 0 + A this happens. 


V. CONSTRAINTS FROM THE LOW-ACCELERATION 
REGIME 


d 2 x ^ M dx a dx^ 

dr 2 = ~~ 


(64) 


can now be expressed in terms of the coordinate time (i.e. the 
time measured by an observer far from the galaxy): 


d 2 a+ ^ dx a dx@ 1 dx^ t dx a dx@ 

dt 2 a P dt dt + c dt dt dt 


(65) 


In order to determine, at least approximately, the range 
of the combination 0 + A for which the PN expansion re¬ 
mains perturbative and the agreement with observations of 
galactic rotation curves (as well as with the observed Tully- 
Fisher relation in disk galaxies and the Faber-Jackson relation 
in ellipticals and clusters) is not ruined, let us consider sys¬ 
tems (galaxies or clusters) with baryonic masses in the range 
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FIG. 2. The bound on /3 + A coming from the requirement that the theory reproduce the rotation curves of galaxies, combined with constraints 
from the high-acceleration regime (solar-system tests, absence of Cerenkov vacuum radiation, BBN, pulsars, classical and quantum stability). 


M = 10 10 — 10 14 Mq. Note that for these masses, the ra¬ 
dius ro marking the onset of the MOND effects lies well out¬ 
side the half-light radius given by Eq. m so our calcula¬ 
tions (which assume that r 0 is larger than the size of the lu¬ 
minous component of the system) do hold, at least to first ap¬ 
proximation. One cmcial ingredient to calculate the impact of 
the 1PN terms on the rotation curves is, as can be seen from 
Eq. ( [69] ), the accretion rate of IGM gas onto the galaxy. A 
very rough estimate for this quantity is M ~ M/tu, where 
i/r ~ 1.4 x 10 10 yr is the Hubble time. 

A useful measure of the impact of the 1PN terms on the 
rotation curves is given by the fractional deviation 

e{M,f,P + A) = |^ >1PN /^ )N — l| ■ (70) 

Clearly, this quantity is a function of /3 + A, but also of the 
galaxy’s mass M and the orbital radius r. Since what is of in¬ 
terest to us is the range of j3 + A for which e(M, r, /? + A) is not 
“too large”, we can marginalize over M £ [10 10 : 1O 14 ]M 0 , 
and over r. For the latter, we marginalize over a range 
spanning from r = ro (M) (the distance from the center at 
which MOND effects become dominant) up to the virial ra¬ 
dius r = r v ; r (M) (at z = 0) of the ACDM halo corresponding 
to the galaxy under consideration. This choice is justified be¬ 
cause rotation curves are measured well beyond the galaxy’s 
luminous part, deep into what in the ACDM model is the dark- 
matter halo region. In order to estimate r v - lr (at z = 0), we 
use J76MI 

( M \ 1/3 

r v ir = -TT- iril3 , /f Mpc, (71) 

V/b x 5.5 x 1O 13 M 0 / 

where /t, « 0.17 is the baryon fraction in the ACDM model. 

In order to identify the range of fj + A for which 1PN terms 
“spoil” the agreement with observations, we then consider the 
marginalized fractional deviation 


e(P + X) 


max [e(M, r, /? 4- A)] 


re[r 0 (M) ,r vir (M)] 

Me[io lo -io 14 M 0 ] ’ 


(72) 


and when this quantity exceeds a certain threshold, we con¬ 
clude that the 1PN terms jeopardize the agreement between 
the theory and observations. Assuming a 20% threshold (i.e. 
e = 0.2), we find that in order for the theory to reproduce 
galaxy rotation curves, one must have 

|/3 + A| > 2.5 x 10 -7 . (73) 

This bound is very conservative, e.g. when considering a 30% 
threshold, and marginalizing only over M £ [10 11 : 1O 13 ]M 0 
and r £ [2ro, 0.5r v ; r (M)], a larger region of the (/3, A) plane 
would remain viable, namely \/3 + A| > 2.5 x 10 -9 . Never¬ 
theless, even with the conservative bound given by Eq. m 
a significant region of the (/3, A) plane remains viable when 
one combines that bound with the constraints from the high- 
acceleration regime (c.f. discussion in section|lH]i. This viable 
region of the parameter plane is represented in Fig. [2] 

VI. DISCUSSION: OPEN QUESTIONS AND PROBLEMS 

In this paper, we have introduced a theory that can repro¬ 
duce the MOND phenomenology (and in particular the rota¬ 
tion curves of galaxies) at low accelerations (i.e. low ener¬ 
gies), and which reduces to khronometric theory/Horava grav¬ 
ity at intermediate/high accelerations (i.e. intermediate/high 
energies), thus satisfying experimental requirements such as 
solar-system tests mmm, binary- and isolated-pulsar 
constraints ATI 02), BBN EH 02] M [661, the existence of 
well-behaved black-hole solutions forming from gravitational 
collapse 1431091 . and the absence of gravitational Cerenkov 
radiation l38l . as well as theoretical requirements such as clas¬ 
sical and quantum stability 11331 1401 l65l . 

This transition from a MOND-like theory to khronomet¬ 
ric theory/Horava gravity is achieved by making one of the 
coupling constants of the theory effectively energy-dependent. 
This was first proposed in Ref. ED, but here we generalize 
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that idea by showing that the theory’s 1PN dynamics becomes 
strongly coupled at low accelerations, unless the other (two) 
coupling constants of khronometric theory also have non-zero 
values. In other words, we show that in order to make the 
1PN dynamics perturbative at low energies, the theory can¬ 
not reduce exactly to GR at intermediate/high energies (as 
was conjectured by Ref. |[50l ), but rather to khronometric the- 
ory/Horava gravity. Of course, one clear shortcoming of our 
approach is that it is purely phenomenological at this stage. 
Indeed, we assume that the running of the coupling constants 
is exactly the one that we need to reproduce data/observations. 
It remains to be seen if this running is actually the one pre¬ 
dicted by the renormalization-group flow, but as far as we are 
aware no studies in this direction are available yet. 

Another open question about our approach (and about 
Lorentz-violating gravity in general) is the nature of the mech¬ 
anism preventing the violations of Lorentz symmetry from 
percolating into the matter sector, where they are strongly 
constrained by cosmic-ray/particle-physics experiments. In 
particular, the higher-order operators that are crucial for the 
power-counting renormalizability of Horava gravity must be¬ 
come important at energies < 10 16 GeV to ensure that the the¬ 
ory remains perturbative in the UV. This scale is comparable 
with the energy at which Lorentz violations can be probed in 
the matter sector, thanks to the synchrotron emission from the 
Crab Nebula (54). The percolation of Lorentz violations into 
the matter sector can of course be suppressed at tree level (by 
assuming that matter does not couple directly to the Lorentz - 
violating field), but it naturally reappears due to radiative cor¬ 
rections. To ensure the viability of the theory, a more efficient 
suppression mechanism must therefore be present. Proposals 
include fine-tuning, “gravitational confinement” f58l . “custo¬ 
dial symmetries” (e.g. softly broken supersymmetry I59ll60l ). 
or dynamical emergence of Lorentz symmetry at low ener¬ 
gies in the matter sector, e.g. due to renormalization group 
flows Elia. 

At a more phenomenological level, a pertinent question is 
whether the theory that we propose can explain all cosmolog¬ 
ical data (besides galaxy rotation curves) with no Dark Matter 
at all. This seems unlikely because MOND itself, as men¬ 
tioned earlier, requires some amount of Dark Matter in the 
center of galaxy clusters - with mass roughly twice that of 
observed baryons Q9). As mentioned, however, this “missing 
mass” problem is much less serious than in the ACDM model, 
since one can postulate that this Dark Matter is given by a 
(small) fraction of the “dark missing baryons” predicted by 
BBN and not yet observed. In particular, these dark baryons 
may be in the form of molecular hydrogen l25l . Another pos¬ 
sibility is that the missing mass in clusters may be given by 
neutrinos DU. (Note that the bounds on the neutrino masses 
and families from the CMB do not hold in MOND, at least 
rigorously, as they assume the ACDM model to start with.) 
Also, we recall that without some amount of Dark Matter 
(in baryons or other components), MOND might have a hard 
time reproducing observations of the “Bullet Cluster” l79l . al¬ 
though the interpretation of the data may be more subtle than 
initially thought, since a similar system - the “Train wreck 
Cluster” lt80l - shows a different behavior. 


On scales even larger than those of galaxy clusters (i.e. 
those relevant for type-la supernovae, CMB and large-scale 
galaxy surveys), the full relativistic theory has to be used, 
in order to account for both the background expansion and 
perturbations about it. For a Robertson-Walker (i.e. homo¬ 
geneous and isotropic) background, and assuming that the 
khronon field is aligned with the cosmic time (i.e. that hy¬ 
persurfaces of constant khronon are also ones of constant cos¬ 
mic time), our theory predicts the same Friedmann-Lemaitre- 
Robertson-Walker equations as in GR, with the only differ¬ 
ences being that (i) no Dark Matter is present (except possibly 
the small amount, e.g. in “dark missing baryons”, needed to 
explain galaxy-cluster data); and (ii) the gravitational constant 
differs from the value G'y measured in the solar-system, and 
is given by Gq = Gjv( 1 — a/2)/(l + /?/2 + 3A/2). Given the 
constraints on a, /? and A discussed in this paper, Gc ~ GVy to 
within a few percent. This probably makes it difficult to repro¬ 
duce both type-la and CMB data. Indeed type-la supernova 
observations are only sensitive to the background expansion 
history, and (to first approximation) constrain a linear combi¬ 
nation of the density parameters of matter (fl m ) and cosmo¬ 
logical constant at z = 0. As for the CMB, a detailed 
study of perturbations over the cosmological background is 
needed to predict the details of its angular spectrum, but the 
position of its first peak only depends, to first approximation, 
on the sound speed of the photon-baryon fluid, and on the an¬ 
gular distance to the baryon-photon decoupling. Both these 
quantities are the same in our theory as in the ACDM model. 
Therefore, because the position of the first CMB peak within 
the ACDM model constrains f2 m + Cl a ~ 1 (a constraint al¬ 
most orthogonal to that coming from type-la supernovae), it 
is clear that our model may have a hard time reproducing both 
CMB and type-la supernova data, unless we allow as much 
Dark Matter as in the ACDM model. A more detailed analy¬ 
sis, however, is needed to confirm this, and will be performed 
in future work. Indeed, one may be able to reproduce the data 
without Dark Matter, but by relaxing the assumption that the 
khronon must be aligned with the cosmic time. 

Another possibility comes from the observation that an ef¬ 
fective Dark-Matter component on large cosmological scales 
naturally arises in theories similar to ours, namely in Horava 
gravity with the projectability condition. That is a theory with 
(infrared) action given by Eq. ((5}, but with a = 0 and the extra 
condition (“projectability”) that the lapse N is only a function 
of time (i.e. a M = 0) at the level of the action. More specif¬ 
ically, Ref. I8l1l (c.f. also Ref. 11821 ) showed that such an ef¬ 
fective Dark-Matter component appears in projectable Horava 
gravity if deviations from homogeneity are present on large 
(even super-horizon) scales. It is also well known that the 
solutions to projectable Horava gravity can be obtained from 
solutions to khronometric theory [action given by Eq. <HJ] in 
the limit a —> oo (30l l35l , or equivalently from solutions to 
our theory [action given by Eq. ( fTIT )] for x ocpl This can 
be shown by following the argument of Ref. |35ll . 


10 We thank Niayesh Afshordi for suggesting this point. 
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Let us then assume that at the scales relevant for galaxies 
and clusters we still have /(a) ss — 2Ao + 2a 2 — 4a 3 /(3ao) 
(so that the results of this paper remain valid), but on larger 
cosmological scales (i.e. even smaller accelerations a —> 0) 
f(a) ~ — 2A 0 + O(ci), so that y oc 1/a diverges as a —>■ 0. 
With this Ansatz, the terms depending on a in the field equa¬ 
tions ([T7}-((26 ]i all vanish when a —> 0, with the excep¬ 
tion of the terms giving the cosmological constant and the 
term Di(xa l ) in the modified Hamiltonian constraint ( fT7| . 
To find the effective Friedmann-Lemaitre-Robertson-Walker 
equations in an inhomogeneous universe, one can take a spa¬ 
tial average of the field equations. In the case of Eq. 
the average of the term Di{xa l ) produces a boundary term 
C, which may not be zero if the universe is inhomogeneous 
on large (even superhorizon) scales. Indeed, this boundary 
term might be interpreted as an effective Dark-Matter com¬ 
ponent, because it has the right scaling with the expansion pa¬ 
rameter A(t), i.e. the spatial average of the modified Hamilto¬ 
nian constraint ( fT7] > yields an effective Friedmann-Lemaitre- 
Robertson-Walker equation A 2 + kc 2 = 8ttGA 2 ( p + pdm)/3, 
with pdm = C/A 3 . Note that this effective Dark-Matter com¬ 
ponent might also improve the agreement of the theory with 
galaxy-cluster observations, which as mentioned above show 
some tension with MOND. 

Finally, another possibility would be to replace the term 9 2 
in the action © with a function of 9 2 . Since 9 is essentially 
given by the Hubble rate for a cosmological background, this 
change may provide enough freedom to fit the background’s 
expansion history, possibly even providing an effective “Dark 
Energy” component. (Note that this is similar to the “general¬ 
ized” Einstein-/Lther theories introduced in Refs. or 

the “K-essence” of Ref. l83l .) Clearly, such a modification 
of the action ( fTTfr may affect the analysis of the 1PN dynam¬ 
ics that we performed in this paper, but the formalism that we 
developed here is readily extensible to that case. 

Of course, all of these possibilities require further detailed 
exploration before one can make any definitive claims about 
them. We will study them, both at the level of the cosmo¬ 
logical background and perturbations about it, in subsequent 
publications. 


APPENDIX: THE PPN PARAMETERS IN THE 
HIGH-ACCELERATION REGIME 


In this Appendix, we show how to solve the 1PN dynam¬ 
ics in the high-acceleration regime, where our theory reduces 
to khronometric theory/Horava gravity. In particular, we con¬ 
firm, as already shown in Ref. [[33, S3, that all the PPN pa¬ 
rameters of khronometric theory are the same as in GR, with 
the exception of the preferred-frame parameters u\ and 0 - 2 . 

At high accelerations (where x = a )> Eq. ( |35| yields the 
usual expression for the Newtonian potential. 


4>n = —Gn J d 3 x' g'j , (Al) 


where we recall that the locally measured gravitational con¬ 
stant Gn is related to the “bare” one appearing in the action 
by Eq. Q. The equations characterizing the 1PN dynamics 
are Eqs. (|42[) and ( |43| , which in the high-acceleration regime 
become 


V ( ^ ^ J = 4t tG nP 

8ttGnPV 2 + 127 tGnP + 2V<f>jv ■ V0 at 

c A L 

2 + /3 + 3A 


V Wi T 2 


2-a 

(£t) 


-(<9 t V 2 w + 8d/(j) N ) 


(A2) 


d,V 2 


16irGpVi 

1—0 


- 2 ( 2+ /- + l 3A )'W~- (A3) 


Before solving them, let us first define the PN potentials Bl: 
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X(x,t) = GN J d 3 x'p(x',t) \x — x! 

h = g„ /dV 


\x — x' 


W, = Gm / d 6 x‘ 


l3 , p{x\ X')i 


\x — X'\ 

$1 = Gn 


\x — X' 


$2 = ~Gn I d 3 x‘ 


$4 = Gn / d 3 x 


3 ,p{x' ,t)<j) N (x' ,t) 


x — x' 

rPjx'A) 
\x — x'\ 


(A4) 

(A5) 
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and recall the following relations among them (4): 


ing of Wi. The solution for got then reads 


V 2 X = —2 <j> N , 

(A10) 

V 2 Vi = -47 TG N pvi, 

(All) 

V 2 $i = -4ttG nP v 2 , 

(A 12) 

V 2 4> 2 = 47 tGnp4>n, 

(A13) 

V 2 $4 = -4nG N p, 

(A 14) 

d,V l = dt<l> N , 

(A15) 

d t V l = -diW\ 

(A 16) 

d t d{K = Wi-Vi. 

(A17) 


Wi diUJ 

~ ^3 + + °( 5 ) 

/3 1 2 + A + 3/3(1 + A) — a(l -)- /3 + 2A) W) 
= 2(/3 — l)(/3 + A) ^ 

a + 5/3 - 3a/3 - /3 2 + A(7 - 2a - 3/3) V, 


+ 


2(/3-l)(/3 + A) 


0(5). 

(A20) 


By using the solution for to and the relations ( |A 1 0[ )-( [A 1 1\ , one 
can then solve Eq. (|A2[) for obtaining 


= <t>N~ 2$i - 2$ 2 - 3$ 4 


(a- 2/3)(2 + /3 + 3A) 2 
2(a — 2)(/3 + A) 4 ’ 


which yields the complete solution for goo at 1PN order: 


Equation (|A3|i can then be written as 


9w = -1-2%-24 

cr c 4 


4 fl +4 + 

,4 ' M 


. <l>, 


C* C 

(a — 2/3) (2 + /3 + 3A) 9 2 X 


(a — 2)(/3 + A) 


0 ( 6 ). 


(A22) 


VV + 2(|^)ftV 2 W = 
- 2(2- a r + ( 2 Y + 3A ]w2x 


(A18) 


Taking the divergence of this equation and using the relations 
above between the PN potentials, we then obtain the solution 
for w, i.e. 


Finally, by performing a gauge transformation t —> t + St 
(with St oc <9 t X), we can write the 1PN metric in the standard 
PN gauge, i.e. 


goo = — 1 — 






+ 6^+0(6) 

(A23) 


goi + a 2 )^| - 

9ij = (l _ 2 “+)^3 + ^( 4 ) 


1 

2 


1 + a 2 



(A24) 

(A25) 


a + /3 + 3A 
2(ff+ A) 3 ' X ' 


(A19) 


which, when replaced back in Eq. (|A181, allows the comput 


where the preferred frame parameters are given, as in 
Refs. |32i40l, by 


ai = 


(A26) 


a 2 — 


4(a - 2/3) 

(3-1 ’ 

(a — 2/3) [—/3 (3 + /3 + 3A) — A + a(l + /3 + 2 A)] 


(a-2)09-1)00 +A) 


(A27) 
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